Intrinsic Half-Metallicity in Modified Graphene Nanoribbons 
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We perform first-principles calculations based on density functional theory to study quasi one- 
dimensional edge-passivated (with hydrogen) zigzag graphene nanoribbons (ZGNRs) of various 
widths with chemical dopants, boron and nitrogen, keeping the whole system isoelectronic. Grad- 
ual increase in doping concentration takes the system finally to zigzag boron nitride nanoribbons 
(ZBNNRs). Our study reveals that, for all doping concentrations the systems stabilize in anti- 
ferromagnetic ground states. Doping concentrations and dopant positions regulate the electronic 
structure of the nanoribbons, exhibiting both semiconducting and half-metallic behaviors as a re- 
sponse to the external electric field. Interestingly, our results show that ZBNNRs with terminating 
polyacene unit exhibit half-metallicity irrespective of the ribbon width as well as applied electric 
field, opening a huge possibility in spintronics device applications. 
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Nanomaterials of carbon, like nanotubes, fuUerenes, 
etc., have been of great interest in condensed- matter and 
material science because of their novel low-dimcnsional 
properties [H, 0]. Over past few decades, cutting edge 
research has been carried out for advanced device inte- 
gration, exploring the electronic and mechanical proper- 
ties of these systems. The recent addition in this journey 
is graphene: a strictly two-dimensional flat monolayer of 
carbon atoms tightly packed into a honeycomb lattice 
[1, 3 . Since its innovation 0, d, 0] , it has made possible 
the understanding of various properties in two dimen- 
sions, by simple experiments and has opened up huge 
possibilities for electronic device fabrications [1, M, [l3l ■ 
A large number of theoretical and experimental groups 
all over the world have gathered on this two dimensional 
platform to search for the "plenty of room" at this re- 
duced dimension 
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Electronic properties of low dimensional materials are 
mainly governed by their size and geometry. Recent 
experimental sophistications permit the preparation of 
finite size quasi one dimensional graphene, named as 
graphene nanoribbons (GNRs) of varying widths, either 
by cutting mechanically exfoliated gra phenes and pat- 
terning by lithographic techniques [12l . 13| or by tun- 



ing the epitaxial growth of graphenes 

Differ- 
ent geometrical terminations of the graphene monolayer 
give rise to two different edge geometries of largely vary- 
ing electronic properties, namely, zigzag and armchair 
graphene. Several theoretical models, e.g., tight-binding 
model within Schrodinger 16 JTtI . [l8l|. Dirac formalism 
for mass less fermions [l3,[Ml2lf7 density functional the- 
ory (DFT) etc. have been applied to explore the elec- 
tronic and band structure properties of GNRs. There 
exists a few many-body studies, exploring the electronic 
and magnetic properties of these systems . 

DFT studies suggest that, the anti-ferromagnetic quasi 
one-dimensional (ID) zigzag edge graphene nanoribbons 
(ZGNRs) show half-metallicity at a finite external elec- 




FIG. 1: (Color online) Schematic representation of the sys- 
tems considered. The upper device is fabricated with hydro- 
gen passivated 8-ZGNR (8 zigzag chains) without any doping. 
The lower device depicts the doping positions in 8-ZGNR, 
replacing two middle zigzag carbon chains by isoelectronic 
boron-nitrogen chains. Doping concentration increases gradu- 
ally along the directions of solid arrows. Electric field (dashed 
arrows) has been applied along the cross-ribbon width, which 
is the y-axis direction. 



trie field across the ribbon width within both local den- 
sity approximation (LDA) 



24| and generalized gradient 



approximation (GGA) [25| . Half-metallic materials show 
zero band gap for electrons with one spin orientation and 
insulating or semiconducting band gap for the other, re- 
sulting in completely spin polarized current. This in- 
teresting property has been observed in some materials 
like Hcusler compound s 1201 , manganese perovskites 27| , 
mctal-DNA complexes [28[ etc. Modification of edges and 
systematic doping indeed change the electronic proper- 
ties of GNRs and thus give new inroads in fine tuning 
of its band gap, controling the nature of the spin polar- 
Experimental progresses have 
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ization [29|, | 

enabled realization of such materials to a great extent 
[H, [13, H [H, Hi. In this letter, we focus on the 
effect of chemical doping on edge passivated ZGNRs. 
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FIG. 2: (Color online) The Density of states as a function 
of energy, scaled with respect to the Fermi energy for the 8- 
ZGNR with (al) 2, (bl) 4 and (cl) 6 zigzag carbon chains 
replaced by the zigzag boron-nitrogen chains at the middle of 
the nanoribbon. The spin density profile of these systems are 
shown in (a2), (b2) and (c2) respectively. 



By chemical doping, wc mean the replacement of car- 
bon atoms by some other elements, keeping the system 
isoelectronic. 

This can only be achieved if two carbon (ls^2s^2p^) 
atoms in ZGNRs are replaced by one nitrogen 
(ls^2s^2p'^) and one boron (ls^2s^2p^) atom. The effect 
of these chemical dopants has already been studied both 
experimentally and theoretically in case of nanotubes 
37|. It has been observed that, the ratio of carbon atoms 
to these dopants has significant impact on the magnetic 
behavior of the nanotubes and increse in boron and ni- 
trogen concentration gradually opens up the gap between 
the valence and conduction band. Here we start with 
undoped hydrogen passivated 8-ZGNR (8 zigzag chains) 
and incorporate the dopants by replacing the two middle 
zigzag carbon chains by two zigzag boron-nitrogen (B- 
N) chains in such a way that, each C-C unit is replaced 
by one B-N unit, resulting in a situation where the mid- 
dle of the ZGNR looks like as a polyborazinc (Figure 1). 
We then systematically increase the doping concentration 
by replacing the two adjacent zigzag carbon chains, just 
next to the polyborazine at the middle and keep doing 
so upto the edges. Ultimately we end up at the zigzag 
boron nitrogen nanoribbon (ZBNNR), where all the C-C 
units arc replaced by B-N pairs. Although the ZGNRs 
and ZBNNRs are isoelectronic, they exhibit completely 
different electronic properties. As like bulk hexagonal 
boron nitride, quasi one-dimensional ZBNNRs show wide 
gap insulating behavior 
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|. In this letter, we 
consider chemically doped ZGNRs with varying widths 
and study their electronic properties with and without 
external electric field. 

We perform first principle periodic calculations based 
on density functional theory adopting the SIESTA pack- 



age [411 ]. Within generalized gradient approximation 
(GGA), considering Perdew-Burke-Ernzerhof (PBE) ex- 
change and correlation functionals 42] with a double zeta 
polarized (DZP) basis set, we have performed spin polar- 
ized calculations. The pseudopotentials are constructed 
from 3, 4 and 5 valence electrons for the B, C and N 
atoms, respectively. GGA approximation takes into ac- 
count the semilocal exchange correlations which have 
significant impact on low-dimensional spin systems like 
GNRs and BNNRs. We have considered the nanorib- 
bons translated infinitely along the x-axis, as shown in 
Figure 1. We consider the lattice vectors along the y 
and z-axis as 50.0 A and 16.0 A respectively, creating 
sufficient vacuum to avoid the interactions between the 
adjacent unit cells. We have considered a 400 Ry energy 
cutoff for a real space mesh size and a k-point sampling 
of 36 k-points, uniformly distributed along the 1 D Bril- 
louin zone. We have considered nanoribbon widths of 
8, 12 and 24 zigzag chains with varying concentration 
of doping and relax all of the structures with both fer- 
romagnetic and anti-ferromagnetic spin orientations. A 
transverse electric field has been applied along the cross- 
ribbon direction, i.e., along the y-axis (see Figure 1) on 
the ground state structures to mimic the band gap engi- 
neering for both the spins. 

All the systems show zero ground state magnetic mo- 
ment. The band structure analysis shows that, the un- 
doped ZGNRs are narrow gap semiconductors with simi- 
lar gaps for up and down spins. On application of ex- 
ternal electric field, one of the spin gaps closes at a 
certain critical field strength, making the system half- 
metallic. The direction of applied field governs the na- 
ture of the conducting spin channel 2J, |25[ . On the other 
hand, the ZBNNRs are known as wide gap insulators 
and the gap closes on application of electric field with- 
out any spin polarization of the current [s^ . These are 
two extreme cases. Gradual replacement of zigzag car- 
bon chains at the middle of the ZGNRs with zigzag B-N 
chains takes the ZGNRs to ZBNNRs and the electronic 
and band structures change accordingly. In Figure 2, we 
have shown the density of states (DOS) of 8-ZGNR with 
replacement of (al) 2, (bl) 4 and (cl) 6 middle zigzag 
carbon chains. The undoped 8-ZGNR shows band gap of 
0.5 eV for both the spins [l^. Replacement of 2 middle 
zigzag carbon chains gives non-zero DOS near the Fermi 
energy for one of the spins, reducing the gap significantly 
but keeping the other spin channel semiconducting. The 
system with 4 zigzag B-N chains at the middle makes one 
of the spin channels conducting. But there is an increase 
in the gap for the other spin channel. Interestingly, the 
system with 6 zigzag B-N chains at the middle and two 
zigzag carbon chains at the edges shows decrease in both 
the spin gaps, a distinct behavior different from the other 
two cases. The DOS analysis in first two cases suggest 
that the increase in doping concentration closes one of 
the spin gaps and opens up for the other spin channel. 
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FIG. 3: (Color online) The density of states as a function of 
energy, scaled with respect to the Fermi energy for (a) 8, (b) 
12 and (c) 24 ZBNNRs with terminating polyacene unit. 



These observations clearly indicate a deviation from the 
symmetric spin distribution in bipartite ZGNRs. 

The ground state of ZGNRs with zero magnetic mo- 
ment, show antiparallel spin alignment between nearest 
neighbours, a typical nature of bipartite lattice which has 
earlier been noted [IS, 0, IH] . Although all the atoms in 
one edge are aligned ferromagnetically, they show anti- 
ferromagnetic coupling within two edges. However, in the 
ZBNNRs, both the spins prefer to mix homogeneously 
throughout the whole lattice without any dominance of 
a particular spin on any site. Figure 2 shows the spin dis- 
tributions of the intermediate systems, i.e., zigzag boron- 
nitrogen chain doped ZGNRs. The figure suggests that, 
whenever carbon atoms find nitrogen or boron as near- 
est neighbor, the bipartite symmetry breaks. Moreover, 
due to the Lewis acid character, the boron atoms pull 
the electrons from adjacent carbon atoms, resulting in 
a charge transfer from carbon to boron and creating a 
potential gradient across the ribbon width (see Figure 
2). Boron and nitrogen face different spins on the ad- 
jacent carbon atoms in either side of the doped zigzag 
lines. Thus unsymmetrical spin diffusion on either side 
destroys the spin symmetry, which is reflected in the 
DOS calculations. The spin distribution of the system 
with 6 zigzag B-N chains at the middle and terminating 
polyacetylene chains (Figure 2 (c2)) shows localization of 
spins on one of the edges only, resulting in a completely 
inhomogeneous electronic structure resembling neither 
ZGNRs nor ZBNNRs. Note that, ah the other doped 
systems are mixtures of polyborazine and polyacene units 
which are the constituents of ZGNRs and ZBNNRs re- 
spectively and thus exhibits mixed properties. Infact, 
carbon dopants in this system (ZBNNR with terminat- 
ing polyacetylene chains) insert finite density of states for 
both the spins within the large insulating gap of ZBNNR. 
To estimate the contribution of the terminating carbon 
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FIG. 4: (Color online) Contour plot of the density of states as 
function of energy, scaled with respect to the Fermi energy for 
(a) conducting and (b) insulating spin channels of 8-ZBNNR 
with terminating polyacene unit with varying electric field. 



atoms, we perform calculations for projected density of 
states (pDOS), which indeed shows that, the DOS near 
Fermi energy come solely from the edge carbon atoms. 

pDOS analysis for all the systems show that the maxi- 
mum contribution comes from the carbon atoms and the 
extent of their contributions reduce from the edges to- 
wards the middle. However, the carbons, attached to the 
boron atoms and subsequent carbon atoms on the same 
side of doped ZGNRs shows higher contributions than 
that of the other side, which is attached to the nitro- 
gen atoms. This observation unambiguously proves that, 
the half-metallicity in boron-nitrogen doped ZGNRs orig- 
inates from the charge transfer from the carbon to boron 
atoms. This can be further explained from the bare on- 
site Coulomb repulsion values of carbon (~ 9.66 cV), 
boron (~ 8.02 eV) and nitrogen ('^ 14.46 cV), which has 
been estimated as the difference between the ionization 
potential and electron affinity of the bare atoms. Due to 
reduced onsite Coulombic replulsion, the electrons pre- 
fer to flow from carbon to boron, whereas the charge 
transfer from carbon to nitrogen is hindered by the large 
Coulombic repulsion on nitrogen sites. Although these 
atomic values renormalize in case of complex systems like 
nanoribbons, their relative differences remain almost un- 
affected. 

Calculations with varying ZGNR width with all possi- 
ble doping concentrations result in similar observations. 
Our study reveals that the systems with two zigzag car- 
bon chains on either edges and all zigzag B-N chains at 
the middle shows half-metallicity irrespective of the rib- 
bon width, which has been shown in Figure 3. pDOS 
analysis for all the widths confirms the Lewis acid char- 
acter of boron for the observed half-metallicity. The in- 
sulating gap for the other spin channel decreases with 
an increase in ribbon width, manifesting the characteris- 
tics of both ZGNRs and ZBNNRs in absence of doping 



39j. In Figure 4, we have shown the variation of 
the band gap for both the spin channels for 8-ZBNNR 
system with terminal polyacene units. Application of 
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FIG. 5: (Color online) The spin polarized band structure (5 
valence + 5 conduction) of 8-ZBNNR with terminating poly- 
acene unit, as obtained from (a) SIESTA and (b) VASP cal- 
culations for (i) up and (ii) down spins separately. The Fermi 
energy (Ep) is set to zero. 



electric field on ZBNNRs with polyacene edges for all 
widths shows that, the half- metallic behavior sustains 
over even sufficiently large e-field strength. The mag- 
nitude of the semiconducting spin gap (> 0.3 eV) sug- 
gests the possibilities of experimental realization of such 
half-metallic anti-ferromagnets even at room tempera- 
ture (0.025 eV). 

We have also verified that, the intrinsic half-metallicity 
as obtained for ZBNNRs with terminating polyacene unit 
from SIESTA with localized basis remains qualitatively 
same with extended plane augmented wave 4^ basis 
within DFT package VASP [43. In Figure 5, the spin 
polarized band structure results for 8-ZBNNR from both 
SIESTA and VASP at the same level of calculations are 
shown. As can be seen, although the gap for the insu- 
lating spin channel decreases for extended basis states, 
it remains far above the room temperature, highlighting 
our claim for intrinsic half-metallicity in ZBNNRs. 

In conclusion, we have studied the ZGNRs with sys- 
tematic chemical doping, replacing the zigzag carbon 
chains with zigzag boron-nitrogen chains in the middle of 
the nanoribbons, keeping the systems isoclcctronic with 
the undoped ZGNRs to search for better half-metallic 
materials. Our study suggests that, the nanoribbons 
with terminating polyacene unit and all zigzag boron- 
nitrogen chains at the middle has the potential to act 
as a half-metallic anti-ferromagnet for all widths. We 
have identified the Lewis acid character of boron, that 
results in charge transfer form carbon to boron to invoke 
half-metallicity in graphenc based materials. This prop- 
erty sustains over a large external electric field strength 
even at room temperature, opening huge possibilities for 
the realization of graphenc based spintronic devices with 
complete control over the spins. 
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